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Abstract 

Symmetry and reciprocity constraints on polarization state of the field diffracted by 
gratings of quasi-planar particles are considered. It is shown that the optical activity 
effects observed recently in arrays of quasi-planar plasmonic particles on a dielectric 
substrate are due to the reflection of the field at the air-dielectric slab interface and are 
proportional to this reflection coefficient. 

1 Introduction 

Two-dimensional periodical grids of thin planar metal particles of various complex shapes 
show interesting electromagnetic properties and find applications in radio and microwave 
engineering as frequency selective surfaces and polarization transformers (e.g., [1]). With 
recent advances in nanofabrication, optical properties of such structures start attracting 
considerable interest of researchers. In particular interest are 2D-chiral inclusions, that do 
not possess mirror symmetry under two-dimensional reflections with respect to a line in 
the particle plane. Various gammadion shapes have been recently studied theoretically and 
experimentally in the optical region, e.g., [2, 3, 4]. 

The symmetry of these shapes has an interesting analogy with magnetized ferrites and 
with three-dimensional (volumetric) chiral particles. The gammadion (swastika) geometry 
defines an axial unit vector normal to the particle plane, whose direction is determined by 
the tilt direction of the particle arms and the right-hand rule. Similarly, in the case of a 
magnetized ferrite sample, there is a selected direction defined by the bias magnetic field. 
At first sight, this suggests that in arrays of 2D-chiral particles there can be electromagnetic 
phenomena analogous to the Faraday rotation of the polarization plane. And indeed, polar- 
ization rotation (optical activity) has been observed in arrays of swastika- shaped particles on 
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a dielectric [2, 3, 4]. However, in [3] this phenomenon was mistakenly attributed to broken 
time-reversal symmetry or nonreciprocity (note that the Faraday effect is nonreciprocal). 
Later it was experimentally shown that the rotation of the polarization plane of transmitted 
waves was reciprocal [4]. Moreover, in [4] the dielectric substrate was found to be the key 
factor for the optical activity in quasi-planar arrays of gammadions because adding a sub- 
strate would effectively turn an array of 2D-chiral objects into a 3D-chiral structure which 
would exhibit the usual reciprocal optical activity. 

Despite apparent similarity, polarization rotation phenomena arc fundamentally different 
in reciprocal and nonreciprocal systems. In the case of magnetized ferrites, time reversal 
operation reverses the direction of the axial vector defining the bias field direction, while in 
the case of arrays of reciprocal particles, time reversal does not affect the direction of the 
axial vector showing the particle "rotation sense". Any system formed by metal and dielectric 
inclusions of any shapes is reciprocal, because there is no time-odd physical parameter there 
and the Maxwell equations themselves are invariant under time reversal operation. 

It is well known that polarization plane rotation exists also in reciprocal media. Optical 
activity is the main manifestation of three-dimensional chirality of molecules or artificial 
inclusions forming chiral and, more generally, bi-anisotropic materials [5, 6]. This effect 
is governed by the chirality parameter (the trace of the magnctoclcctric coupling dyadic), 
which vanishes for planar particles of arbitrary shapes, because any 2D shape can always 
be superimposed with its mirror image if rotations in 3D space are allowed. In fact, the 
magnetoelectric coupling dyadic of planar particles is antisymmetric, and composites formed 
by such particles belong to the class of omega media [6]. This apparently imposes restrictions 
on optical activity effects in arrays of 2D particles, even if their shape is chiral in the two- 
dimensional space. 

Electromagnetic properties of 2D-chiral arrays in the microwave region were studied ear- 
her (e.g. [7, 8]) and these structures were even patented [9]. However, to the best of our 
knowledge, no systematic consideration of fundamental limitations on electromagnetic scat- 
tering from general grids of quasi-planar particles has been made. Important results were 
obtained in [7, 10, 11] from the symmetry considerations, but those results concern restric- 
tions on the form of the constitutive tensors of composite media but not the diffraction 
phenomena. The theory of [12] concerns only arrays in free space and does not include the 
effects of the substrate. 

In this paper we make a systematic study of reciprocity and symmetry restrictions on 
plane-wave diffraction by 2D arrays of arbitrary shaped quasi-planar particles. By quasi- 
planar particles we understand particles of a finite thickness but with uniform properties 
over the full particle thickness. The arrays can be positioned on planar dielectric substrates, 
and we will see that the presence of the substrate is of key importance for optical activity 
phenomena in these structures. 

2 Problem formulation 

We consider an infinite two-dimensional array of arbitrary shaped (in the array plane) parti- 
cles backed with a slab of a homogeneous dielectric. The thickness of the particles is assumed 
to be finite, but the structure is uniform over its thickness. This implies that the array has 
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a plane of symmetry that itself can be regarded as the array plane. In practice, the particles 
are usually made of metal, but the theory which we develop here applies to any reciprocal 
isotropic material. In case of metal inclusions, the general theory developed here applies 
both to microwave frequencies and to the visible range, because particular electromagnetic 
parameters of metal particles do not affect restrictions which follow from symmetry and reci- 
procity considerations. We assume that the array particles are not touching the dielectric 
substrate, while the separation of the dielectric and the particles can be arbitrary small. This 
assumption allows us to split the entire structure in two parts: the grating and the slab and 
consider them separately. An example of the problem geometry is illustrated in Fig. 1. 



dielectric layer 




Figure 1: An array of metal (plasmonic) particles backed by a dielectric slab, n is the unit 
vector normal to the array plane. 

The array is illuminated by a plane wave of arbitrary polarization and arbitrary propa- 
gation direction^ incoming from the half space z < 0. We are interested in the polarization 
states of the transmitted and reflected waves^ and how it relates to the polarization of the in- 
cident wave. We will attack this problem with a rather general approach taking into account 
symmetry and reciprocity restrictions. 

3 Symmetry and reciprocity restrictions 

Let us start from considering the array of particles separately from the dielectric substrate. 
Under our assumptions there exist two planes in space which enclose the grating entirely 

""^ Evanescent- wave excitation is not excluded from this analysis. 

^The case of multiple propagating modes created by the array is not excluded: In that case the following 
restrictions apply to every excited mode. 
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Figure 2: Definitions of the incident and scattered waves for a quasi- planar grid. The arrows 
show the normal to the interfaces components of the wave vector. 



in between themselves. We split the tangential electric fields on these two planes into the 
incident fields exciting the grating and the fields scattered by the grating (see Fig. 2) and 
introduce the following scattering matrix: 

fErW] (Su Su] (Et{r)] 
lE+(r)J [s,^ S,,} ' [E,-(r)J " 

Here Ej|'"(r) and E2 (r) denote the vectorial complex amplitudes of the tangential components 
of the incident electric field on the first and the second surface, respectively. E]"(r) and 
E2 (r) are the respective scattered fields, and Smn are scattering operators that connect the 
tangential components of the scattered and incident electric fields at surfaces 1 and 2. 
The fields at the two surfaces can be expanded into plane waves 

E(kt) = JJ E(r) e-'^' ''dh, (2) 

where the integration is taken over the surface plane, E can be any of the fields defined above, 
and kt is the wave vector component tangential to the planes. We use the time dependence 
of the form exp(—iu}t) conventionally used in optics. 

This expansion gives us a representation of the operators Smn in the basis of plane waves. 
The corresponding matrix elements of the operators Smn in this basis are dyadics Smn(p, q): 

lmn(p, q) = // e"^^" Smn e'^'^dFv. (3) 

If the grating has a center of symmetry and the planes where we define the fields are located at 
symmetric positions with respect to the same center, the dyadics s^n(P) ^} obey the following 
restrictions: 

^11 (P, q) = ^22(-p, -q), ti2(p, q) = i2i(-p, -q). (4) 

Indeed, applying the coordinate transform r — r where the origin is at the symmetry 
center, we come to the same structure with the first and the second interfaces interchanged. 
On the other hand, the components of polar vectors change signs under this transformation, 
and we obtain (4) . 
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If the grating plane is itself a plane of symmetry of the system, then, more obviously, 

Sii = S22, S12 = S21, (5) 

or, which is the same, 

^11 (p, q) = ^22 (p, q), ^12 (p, q) = ^21 (p, q)- (6) 

To get this result, one has to apply a transform that changes the sign of the radius vector 
component normal to the grating plane. Such a transformation interchanges the interfaces 
as above, but preserves the signs of the tangential components of polar vectors. 

Furthermore, if the grating has rotational symmetry of order m > 2 with respect to an 
axis parallel to n, the following additional relation holds: 

^ijil • p, J • q) = J • (p, q) • J ■ (7) 

where ^ = exp[(27r/m)[n x I^]] is the dyadic of rotation in the grating plane. It is the unit 
dyadic in the same plane. 

In a reciprocal system the reciprocity theorem holds. To apply this theorem we assume 
that the external fields are created by sheets of surface electric currents Ji,2 placed at the 
first and the second interfaces, respectively. Then, at the interfaces 

Et(v) = -I^-^, E,-(r) = -^^i:i2W, (8) 

where Z^, is the free-space wave impedance operator. In the plane-wave basis it is 

— — — A ~ qq/^^ 

^w(p, q) = (27r)2 S{q - p) ^„(q), ^„(q) = rj \ (9) 

a/1 - q^lk^ 



where r\= y/io/^o, ^ = cuy^^o/Io- The operator is diagonal and symmetric as is directly 
seen from (9): 

^w(p, q) = ^w(-q> -P); ^w(p, q) = 0, if p 7^ q. (10) 

Because the external field itself satisfies the reciprocity theorem, we can write for the scattered 
fields 

E--Jid'r+ yy"E+-J2d'r = • Jr d'r + E+ • d'r. (11) 

where Ji,i', J2,2' are the external currents and ^x,v^ -^2^2' scattered fields in a pair of 

separate excitation scenaria. 

The currents Ji^2, Ji',2' in (H) can be arbitrary chosen and the corresponding electric 
fields can be expressed from (1) and (8). If, for example, J2 = J2' = 0, we get from (1) 
and (8) 

Ji • (^11 • ^w) • Ji' rf'r = Ji, • (^n • ^w) • Ji dS, (12) 
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which in the operator sense means that transpose does not change the operator in braces. 
Therefore, considering all possible cases, we get for a reciprocal grating: 

-S*!! • — ■ Sji, S22 ■ Zyf — Zyj ■ Sj2-, Si2 ■ Zyf — Z^ " 5^^. (13) 

Here we use the fact that Z^ — Z^ as given by (10). Relations (13) can be also understood 
as transpose rules 

Sji — Z^ • Sn ■ Z^, SJ2 — Z^^ • S22 ■ Zyf, Sj^ — Z^^ ■ S12 • Zy,. (14) 
Eq. (14) can be rewritten in terms of the matrix elements of the corresponding operators: 

^n(-q> -p) = ^w^(p) • ^ii(P' q) • ^w(q), 

S(-q, -p) = ^w^(p) • ^22(p, q) • ^w(q), (15) 

^2i(-q, -p) = ^w^p) ■ ^i2(p, q) ■ ^w(q). 

From (15) one can see that for the waves propagating normally to the grating (when p = 
q = 0) reciprocity simply requires 

^n(0,0)=^n(0,0), ^22(0,0)^^22(0,0), ^12(0, 0) = ^21(0, 0). (16) 

Let us note that when compared to the standard mode coupling parameters in the waveguide 
theory (S-parameters), our scattering dyadics SmniP, q) differ in that sense that, first, we do 
not decompose the fields into TE and TM modes, and, second, we do not separately normalize 
the modes, instead we deal with a simple plane-wave basis applicable to all types of modes. 
This allows for a dyadic formalism but results in additional impedance terms in (15). 



4 Multiple layers 

When there are several layers stacked one on top of another (in our case there are two 
layers: the grating and the dielectric) the direct application of (S-matrix (1) becomes tedious. 
Instead, it is better to apply an approach based on transfer matrices. A transfer matrix 
connects pairs of fields given at separate interfaces: 



(Kir)] 
lEr(r)J 



Til 
T21 



T12 
T22 



E2-(r) J 



(17) 



The total transfer matrix of a layered structure is a product of transfer matrices of the layers. 
The operators Tmn can be expressed trough Smn and vice versa: 



Tn 

Sn 



n-l 
'^21 


^12 


— "'^'21 


T21- 


f-1 

^11 ! 


'S'12 = 



-'22, 



^21 — Sii ■ S21, T22 — S12 — Sii ■ S21 ■ S22, (18) 



T, 



22 



T21 ■ T-^i^ ■ T12, 1521 — T-^i^, S22 — ^T^i^ ■ T12. (19) 

From these rather complicated relations one can see that the symmetries considered in the 
previous section generally do not result in simple constraints on (when compared to 
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Smn)- This explains why we have started our formulation from the matrix of Smn operators. 

One exception is the rotational symmetry rule (7) which also holds for T„j„. 

However, there is one important case which results in simple constraints on the elements 
of the transfer matrix. It is the case when the structure is reciprocal and has a plane of 
symmetry, i.e. its Smn operators satisfy both (5) and (14). In this case after some operator 
algebra we obtain 

^11 = ^w"^ ■ ^11 ■ -^w, = Z^^ ■ T22 ■ Z^, = ~-^w^ ■ T21 ■ Z^. (20) 

In the opposite direction, if the operators Tynn satisfy (20) then the operators Smn are such 
that 

S\\ = Z^ ■ S22 ■ Z^, = Z^ ■ S12 ■ Z„, S21 = Z^ ■ S21 ■ Z^. (21) 

Eqs. (20), (21) can be rewritten in terms of the matrix elements of the corresponding operators 
analogously to (15). 



5 Applications 

In the previous sections we formulated a number of symmetry and reciprocity constraints on 
scattering and transfer matrices. Let us consider some consequences of these constraints. 



5.1 Restrictions on optical activity 

First, it is easy to see that the combination of the reciprocity and the central symmetry in 
a structure forbids any optical activity effects for normally incident plane waves (both in 
transmission and reflection into the same normally propagating mode). Indeed, from (4) 
and (16) we get 



Sll (0,0) =522(0,0) 
^21(0,0) =^12(0,0) 



^i(0,0) 

f 1(0,0) 



^2(0,0), 
^2(0,0), 



(22) 



i.e. both transmission and reflection dyadics are fully symmetric. For optical activity to occur 
one needs antisymmetric terms. We have the same result if the central symmetry is replaced 
by the planar symmetry (6). Moreover, if the structure also has rotational symmetry (7) 

with m > 3 then s^n(0, 0) = s^^/t, i.e. the structure does not change the polarization state 
of the transmitted and reflected fields at all. 

Next, let us look at the original problem of an array of quasi-planar 2D-chiral particles on 
top of a dielectric slab. To obtain the total transfer matrix of this system one has to multiply 
the transfer matrix of the grating by the transfer matrix of the slab: 



I rp tot I 

I mn J 





fl2 


( 


1^21 


T22 ^ 





r- 



(23) 



where T and V. denote the transmission and refiection operators of the slab, respectively. 
The matrix elements of these operators are dyadics t(p, q) and p(p, q): 



t(p, q) = (27r)2 5(q - p) T(q), p(p, q) = {^^n f 5(q - p) p(q). 



(24) 
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where T(q) and p(q) are dyadic transmission and reflection coefficients for a dielectric slab 
as functions of the transversal wave vector of the incident wave (k^ = q). The explicit 
expressions for them can be found in textbooks on electromagnetics. 

We are interested in T/^"' because it is connected with the transmission operator of the 
total system: ^2*1* = {f-^T)~^- I^om (23) we get 

= fn • f-^ + fi2-n- f-\ (25) 

or, rewriting it in terms of the matrix elements, 

— tot H — 

hi (P, q) = ^11 (P, q) • ^ (q) + h2 (p, q) • p(q) • r (q) . (26) 

Let us now look at the transpose of the operator 7\\°* in the case when the grating is reciprocal 
and has a plane of symmetry. In this case 

tot)T ^ f-1 . ^-1 . f^^ .z^-f-^.-R. z-^ . . (27) 

where we used (20) and the fact that the transmission and reflection operators of a dielectric 
slab are diagonal and symmetric: T = T^, TZ = TZ^ . In terms of the matrix elements (27) 
becomes 

= totT ^ r — ^ — 1 

^11 (-q, -P) (p) • (p) • iii(p, q) - p(p) • (p) • i2i(p, q)J • ^w(q)- (28) 

Eqs. (26) and (28) greatly simplify for the plane waves propagating along the normal: 

=tot -, = — 

tn(0,0) = T-^in(0,0)+pr-iti2(0,0), (29) 

=totT - - 

^11 (0,0) = T-Uii(0,0)-pr-^t2i(0,0), (30) 

where p and r are scalar reflection and transmission coefficients for a plane wave normally 
incident on a dielectric slab. From the last relations one can obtain an explicit formula for 

= tot 

the antisymmetric part oit^-^ (0, 0): 

= tot =totT - - = =T 

hiM-hi (0,0) _ ti2(o,o) + t2i(o,o) _ ti2(o,o)-ti2(o,o) 

2 2^ 2^ ■ ^^^^ 

From here one can see that adding a dielectric slab to a reciprocal symmetric grating can 
indeed result in a system with a non-symmetric transmission dyadic, i.e. one can realize 
an optically active system this way. The higher is the reflection, the higher is the obtained 
optical activity. The physical reason for this is that the dielectric layer positioned on one 
side of the array breaks the mirror symmetry of the array, and the whole structure becomes 
chiral in the three-dimensional space. 

Furthermore, we see that the optical activity of this kind cannot be achieved with gratings 
that have ^12(0, 0) -|-t2i(0, 0) = 0. For example, this holds for any planar (thickness tends to 
zero) grating made of a non-magnetic material. Indeed, for such a grating 

5*22 = 5*11, 5*21 = Si2 = 1 + Sii^ (32) 
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therefore 

fi2 = -(1 + ^n)-^ ■ = -(1 + ^n')-' - -^ii • (1 + Sn)-' = -f^i. (33) 

This observation leads us to an additional conclusion that the obtained optical activity de- 
pends also on the thickness of the quasi-planar grating. Increasing the thickness (up to a 
certain limit related to the wavelength) one can increase the term ti2(0, 0) +t2i(0, 0) in (31). 

5.2 An array of swastika-shaped particles 

As an illustrative example let us consider a square array of swastika- shaped metal particles. 
It was found experimentally that such a grating when backed with a dielectric slab rotates 
the polarization of normally incident plane waves [4] . When there is no substrate, there will 
be no rotation. Although this follows from our general considerations given above, let us 
show the same using simpler arguments. 




Figure 3: A swastika-shaped particle excited by an external electric field E. Here p is the 
induced electric dipole moment. 

Consider a swastika- shaped metallic particle depicted in Fig. 3. The particle has rotational 
symmetry or order 4, it has a center of symmetry, and it is reciprocal. Suppose that the 
particle is under normal incidence of a linearly polarized plane wave with the electric field 
vector E oriented as shown in the figure. One can decompose the electric field vector of 
the incident wave to two components along x and y axes. Due to rotational symmetry the 
particle reacts equally to both these components. 

The induced electric dipole moment of the particle is a sum of two components p-^ and 
Pj^, and the total dipole moment p = p^, -f- Pj^ of the particle can be written as 

p = ^ • E, (34) 

where a is the dyadic dipole polarizability of the particle. The rotational symmetry requires 
that Qxx — o,yy, ttxy — —ttyx- Ou the other hand, the reciprocity requires axy — ayx, therefore 

(^xy = ciyx = and a — alt, where It is the two-dimensional unit dyadic. This means that 
the induced dipole moment is co-linear with the applied electric field. 

In a square array the interaction of the particles will exhibit the same rotational symmetry, 
therefore for the dipole moment per unit cell of the array and the polarizability in the array 
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we will arrive to the same conclusions as above. The secondary field created by the array in 
far zone is proportional to the dipole moment per unit cell of the array. Its polarization is 
the same as the polarization of the induced dipole moment, and, hence, it is the same as the 
polarization of the incident field. 

We have to note that this simplistic consideration does not include higher multipole mo- 
ments, and, what is more important, it assumes that the effective induced magnetic moment 
per unit cell of the array is either zero (as in swastikas under normal incidence) or does not 
radiate in the far zone (which is the case, for instance, of omega particles under normal inci- 
dence, when the induced magnetic dipole moment is parallel to the propagation direction). 
It is known that one has to provide some means for magnetoelectric interactions to obtain 
optical activity. Adding a dielectric substrate breaks mirror symmetry of a quasi-planar array 
and makes optical activity possible, as we have shown in the previous section. A detailed 
study of how this process can be described in terms of magnetoelectric interactions is outside 
of the scope of this paper. 

6 Conclusions 

Starting from general symmetry and reciprocity (time-reversal invariance) requirements we 
have formulated a number of constraints on the scattering and transmission operators of 
arrays of quasi-planar linear reciprocal particles. We have shown that if a reciprocal array 
is symmetric with respect to its own plane or has a center of symmetry, then this array 
cannot exhibit optical activity for normally incident plane waves both in transmission or 
reflection. For the case of a reciprocal plane-symmetric array backed by a dielectric slab we 
have shown that the optical activity effects for normally incident waves occur due to reflection 
at the air-dielectric interface and the strength of this effect is proportional to the respective 
reflection coefficient. Moreover, even in the presence of a dielectric substrate, to allow for 
optical activity the thickness of the particles in the array should not be neghgibly small. 

The general restrictions on scattering and transmission operators obtained under this 
research apply also to grids that create several propagating modes in reflection and transmis- 
sion (electrically sparse arrays producing diffraction lobes). In the last case, the restrictive 
relations become more complicated and they allow for asymmetries of diffraction phenomena 
under reversal of propagation directions. However, all these phenomena are reciprocal and 
agree with the fundamental time-reversal invariance of the Maxwell equations. 
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